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SUMMARY 


The viscosity and thermal conductivity coefficients are estimated for 
three models of the atmosphere of Jupiter: a Heavy Model consisting of 22% 

helium and 78% hydrogen, a Nominal Model consisting of 11% helium and 89% 
hydrogen, and a Light Model consisting of pure hydrogen. The effect of trace 
elements is neglected. Linearized approximations are used for the transport 
coefficients of the mixtures; these are found to be in almost constant ratio 
to the values for pure hydrogen, independent of temperature. Short basic 
language programs for computing the coefficients are listed. 


INTRODUCTION 


The Jovian atmosphere is presently thought to consist mainly of molecu- 
lar hydrogen, H 2 , perhaps mixed with helium. He, and with traces of methane, 

CH4, and ammonia, NH3. The amount of He seems to be the most uncertain; esti- 
mates range from about 11% by volume for the Nominal Jovian Atmosphere Model, 
to 22% for the Heavy Models, and to 0.0% for the Light Models. The traces of 
methane are thought to be about 0.1% and of ammonia about 0.02%; these are 
significant amounts as far as radiation absorption and emissivity are con- 
cerned , but are trivial amounts as far as thermodynamic and transport proper- 
ties of the atmosphere are concerned. The temperature range of interest is 
from 100 to 1000 K. The thermodynamic properties needed for Jupiter probe < 

design have been previously calculated and tabulated for this range 
(J. Sperans, Ames Research Center, NASA, private communication). Recently, a 
need has appeared for estimates of viscosity and thermal conductivity of the 
Jovian atmosphere, in connection with the subsonic motion of bodies through 
that atmosphere. The present report gives estimates of these quantities based 
on linearized formulas for the transport properties of mixtures and neglecting 
the effect of traces of CH4 and NH3. 


VISCOSITY OF MODEL JUPITER ATMOSPHERES 


Viscosity of Mixtures 

According to rigorous kinetic theory, the coefficient of viscosity of a 
multicomponent gas mixture should be expressed as a ratio of determinants. 


«™«?T er * Hirs ^ hfeld ® r et al. (ref. 1) show that the off-diagonal terms are 

tribitlon P ?n e fh t0 ^ t f 3Ce ® lements of these matrices and that the major con- 
tribution to the viscosity of the mixture can be expressed as 
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The coefficients n lk are defined as 
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where mi is the molecular weight of component i and Q lk i 8 the Chaoman- 
nskog cross sectionfor momentum transfer in collisions between particles i 
and k. The coefficient A ik is the ratio of Q <lf and the cross sections 
appropriate for mass diffusion; it is somewhat greater than unity because 
deeper penetration of the interparticle potential is required ^influence 

?o r 5/3 th^off diaf° r T°T tUm tra r fer ‘ In faCt> if A ik ls assumed equal 
S, f 5“ di 8 1 elements of ^e determinants vanish exactly (ref 1) ♦ 

although this is somewhat too large a value, the effect on Hil frL using ^ 

his approximation is not serious, because terms with m k « mi contribute 
H become ^ equation (2) * With this simplification, the elements 
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and equation (1) can be expressed in the simple linear form 
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if; the effective mean free path A is defined 

Hansen (ref. 2) shows that this value of A* can he cHvon a .<.,,1 . 

leal interpretation as the mean distance traveled by a molecule of type*! 1 * ^ 
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before it loses its excess momentum relative to the surrounding gas medium in 
a series of head-on elastic conditions, since the fraction of excess momentum 
transferred In a single collision between particles i and k is 
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and the fraction 
of f and the 


of excess momentum transferred per unit time is the product 
collision rate 6 ^ summed over all kj 



Multiplying the reciprocal of equation (8) by the mean velocity u of parti- 
cle i leads to the expression of equation (6). Equation (5) can i thus be 
expressed 

n ■ Z ni &°) <9> 
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where the n ± values are the viscosities of the pure components of the mix- 
ture, and AO i s the mean free path in a pure gas of species i at the 
density Thus, 
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Buddenberg and Wilke (ref. 3) empirically find that agreement with exper- 
imental measurements of viscosities of many gas mixtures is improved if each 
term in the summation of equation (10) is multiplied by 0.69; presumably, this 
correction accounts for the deviation of A lk from the 5/3 value assumed. 
Accordingly, the formula used for the calculations of the viscosity of the 
Jupiter atmosphere is 
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* convenient functional ton. u.ed for Q tl 1 „ caking thl. fit 1. 
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Viscosity of Pure Molecular Hydrogen 
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where m is the mass of the H 2 
fer cross section Q for H 2 - Z 


molecule. Accordingly, the 
H 2 collisions is taken to be 


momentum trans- 
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The deviations of equation (14) from the H, 
given by Watson (ref. 4), by Kestin (ref. 5), by 


coefficients of viscosity 
Hilsenrath et al. (ref. 6), 
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and by Keyes (ref. 7) are 
shown In figure 1. The 
parameters of equation (14) 
were chosen to give the 
best fit to Hilsenrath's 
values from 100 to 300 K 
and to Watson's values 
from 300 to 1000 K. The 
latter are expected to be 
the most reliable in that 
range since they have been 
carefully fitted to the 
most recent and most care- 
fully calibrated experimen- 
tal data. Hilsenrath's 
values are based on collision 
integrals calculated using 
a Lennard-Jones 6-12 poten- 
tial, with the force param- 
eters of the potential fitted 
to older viscosity measure- 
ments; these are expected to 
be the best values in the 
low-temperature range. It 
can be seen that the devi- 
ations of equation (14) are 
generally well within 1% of 
the best values. Even the 
old formula given by Keyes 
for the range of 100 to 500 K 
agrees within 3%, 
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Figure 1.- Deviations of equation (14) from 
published viscosities of pure H 2 . 
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and the value document* -»v Kestin at 293 K falls within the 1% range. In 
view of uncertainties in : <e mixture formulas, equation (14) is considered 
adequately precise. 


Viscosity of Pure Helium 

The coefficient of viscosity of pure He is reproduced reasonably well 
over the temperature range of 100 to 1000 K by the formula 
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where ro is the mass of the He atom. Accordingly, the momentum transfer 
cross section Q for He-He collisions is taken to be 
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Figure 2.- Deviations of equation (17) from 
published viscosities of pure He. 


The deviations of equa- 
tion (17) from the He coef- 
ficients of viscosity given by 
Watson (ref. 4), by Kestin 
(ref . 5) , by Kestin and 
Leidenfrost (ref. 8), and by 
Keyes (ref. 7) are shown in 
figure 2. Again the values of 
Watson are considered to be 
the best in the range from 
300 + 1000 K. The formula has 
been fitted to these values 
and to the 100 K value of 
Kestin and Leidenfrost. The 
deviations from Watson's 
values are generally less than 
0.5%; in the low-temperature 
range, where the formula 
bisects the values of Kestin 
and Leidenfrost and of the 
Keyes equation designed for 
the range of 104 to 373 K, 
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deviations are generally less 
than 2%. Again, the uncer- 
tainties are consistent with 
the uncertainties in the mix- 
ture formulas. 


Viscosity of Jovian Atmospheres 

The viscosity of model Jovirn atmosphere mixtures is shown as a function 
of temperature in table 1 for the temperature range from 100 to 1000 K. The 
first column is the temperature; the second column is the viscosity of pure 
H 2 (which corresponds to both the Cool Light Model and Warm Light Model of the 
Jovian Atmosphere) ; the third column is the viscosity of an 89% H 2 - 11% He 
mixture (which corresponds to the Nominal Jovian Atmosphere Model) ; and the 

t h°Jr” vU S 0 ! 1 J y \ 78 , 2 HJ - 22% He ■*««« (which corresponds 

to both the Cool Heavy Model and the Warm Heavy Model of the Jovian Atmo- 
sphere). The Cool and Warm Models have only a different altitude dependence 
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of temperature. The viscosity of the Nominal Model Atmosphere (with 11% He) 

Is almost exactly 25% greater than the viscosity of pure H 2 . The ratio of 
the two is almost Independent of temperature. Similarly, the viscosity of ';he 
Heavy Atmosphere Models (with 22% He) Is almost exactly 47% greater than the 
viscosity of pure 


TABLE 1.- VISCOSITY OP JUPITER MODEL ATMOSPHERES 


Temperature, K 


Light atm. , pP 


Nominal atm.,^ pP Heavy atm. pP 


41.92 

52.98 

62.77 

55.62 

69.92 

82.46 

67.8 

85.08 

100.17 

78.97 

99. C2 

116.52 

89.4 

112.08 

131.87 

99.25 

124.44 

146.42 

108.64 

136.23 

160.32 

117.63 

147.54 

173.68 

126.29 

158.46 

18C.57 

134.66 

169.02 

199.06 

142.78 

179.27 

211.19 

150.67 

189.24 

223.01 

158.36 

198.97 

234.54 1 

165.86 

208.47 

245.81 j 

173.2 

217.77 

256.84 ; 

180.39 

226.88 

267.65 

187.44 

235.81 

278.27 i 

194.35 

244.59 

288.7 ! 

201.15 

253.22 

298.95 







fABLE 2.- VISCOSITY OF JUPITER NOMINAL ATMOSPHERE 


Altitude, km Temperature, K Viscosity, yP Viscosity ratio, 

mlx/H» 


-270.79 

681.77 

-246.93 

637.33 

-224.62 

595.74 

-203.74 

556.83 

-184.23 

520.44 

-165.98 

486.39 

-148.91 

454.55 

-132.96 

424.78 

-118.03 

396.93 

-104.09 

370.89 

“91.06 

346.51 

-78.67 

323.66 

-67.49 

302.25 

-56.86 

282.17 

-46.94 

263.31 

-37.66 

245.59 

-2S.01 

228.92 

-20.97 

213.21 

-13.48 

198.4 

-6.49 

184.43 

0 

171.23 

6.02 

158.77 

11.6 

147 

16.75 

134.3 

21.43 

121.6 

25.79 

116.65 

29.97 

111.7 

34.04 

110.9 

38.09 

110.1 

42.16 

112.35 

46.32 

114.6 

50.57 

117.29 

54.91 

119,98 

59.36 

122.66 

63.9 

125.35 


195.5 

1.256 

186.79 

1.256 

178.45 

1.255 

170.48 

1.255 

162.85 

1.255 

155.56 

1.254 

148.59 

1.254 

141.93 

1.254 

135.55 

1.254 

129.46 

1.254 

123.62 

1.254 

118.03 

1.253 

112.68 

1.253 

107.54 

1.253 

102.6 

1.254 

97.86 

1.254 

93.29 

1.254 

88.88 

1.254 

84.63 

1.255 

80.53 

1.255 

76.55 

1.256 

72.71 

1.256 

68.98 

1.257 

64.85 

1.258 

60.59 

1.26 

58.9 

1.261 

57.17 

1.261 

56.89 

1.262 

56.61 

1.262 

57.4 

1.261 

58.18 

1.261 

59.12 

1.261 

60.04 

1.26 

60.96 

1.26 

61.87 

1.259 


] 
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TAbLE 2.- CONTINUED 


Altitude, km 

Temperature, K 

Viscosity, pP 

Viscosity ratio, 
mlx/Hp 

280.78 

288.16 

295.53 

302.91 

310.29 

200 

200 

200 

200 

200 

85.1 

85.1 

85.1 

85.1 

85.1 

1.255 

1.255 

1.255 

1.255 

1.255 

317.67 

325.05 

332.44 

339.82 

347.21 

200 

200 

200 

200 

200 

85.1 

85.1 

85.1 

85.1 

85.1 

1.255 

1.255 

1.255 

1.255 

1.255 

354.6 

362.39 

370.98 

380.37 

390.57 

200 

221.67 
243.33 
265 

286.67 

85.1 

91.27 

97.24 

103.05 

108.7 

1.255 

1.254 

1.254 

1.254 

1.253 

401.58 

413.39 

426 

439.42 

453.65 

308.33 
330 

351.67 

373.33 
395 

114.21 

119.6 

124.87 

130.03 

135.1 

1.253 

1.253 

1.254 
1.254 
1.254 

468.69 

484.54 

501.19 

518.66 

536.95 

416.67 

438.33 
460 

481.67 

503.33 

140.08 
344.98 

149.8 
154.54 
159.21 

1.254 

1.254 

1.254 

1.254 

1.255 

556.04 

575.96 

596.69 

618.24 

640.6 

525 

546.67 
568.33 
590 

611.67 

163.82 

168.37 

172.85 

177.29 

181.67 

1.255 

1.255 

1.255 

1.255 

1.256 

663.79 
687.8 
712.64 
738.3 

764.79 

533.33 
655 

676.67 

698.33 
720 

185.99 

190.27 

194.51 

198.7 

202.85 

1.256 

1.256 

1.256 

1.256 

1.256 


10 


s? 




TABLE 2.- CONCLUDED 


Altitude, km 

Temperature, K 

Viscosity, pP 

Viscosity ratio, 
mix/H2 

792.11 

741.67 

206.96 

1.257 

820.27 

763.33 

211.02 

1.257 

849.25 

785 

215.06 

1.257 

879.07 

806.67 

219.05 

1.257 

909.73 

828.33 

223.01 

1.257 

941.23 

850 

226.94 

1.258 


THERMAL CONDUCTIVITY OF JUPITER ATMOSPHERES 


Thermal Conductivity of Mixtures 

Theoretical expressions for thermal conductivity of gases and of gas mix- 
tures have not been verified with the same accuracy as for viscosities. This 
occurs partly because a number of intricate corrections are involved in making 
measurements, and also because the heat transfer associated with interna.' 
modes of molecular energy is difficult to treat with the same rigor as that 
associated with the kinetic modes of energy. In the simple Eucken appro: i- 
mation, the internal degrees of freedom are assumed to have a Maxwe’l dipt* 
butlon that is totally independent of the molecule’s velocity. r ' u is : 
tainly a good approximation for the stiff vibrational mod*»f : Az 

molecule at the temperatures of interest, where there little coupling 

between vibrational and kinetic modes during a collisit.. perturbation; 
however, the same may not be exactly true for the rotational modes of the H 2 
molecule. Nevertheless, the Eucken approximation seems to agree just as well 
with experimental data as the more elaborate theoretical models that have been 
devised; consequently, it will be used here in view of its simplicity and the 
fact that the comparist us between experiment and the theoretical model seem 
adequate to the order of accuracy being considered. 

The coefficient of thermal conductivity of a pure gas can be expressed 
with the Eucken approximation 


* ■ I £ ir (I c v + v) < 20 > 

where M is the gram molecular weight, C v is the specific heat per mole of 
the gas associated with its kinetic energy, and C v * is the specific heat 
associated with internal molecular energy. For a Maxwell-Boltzmann distribu- 
tion of velocities, the kinetic specific heat is just 3/2 times the gas 
constant R 

■ ? (2D 
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The internal specific heat is u ^ ^ 

rotational - 1 vibrational contributions; ^ ^ H2 lt 18 the sum of the 


V (He) = 0 


° V * ( " 2) * * [ l + ( ^TT75 t)] 


(22a) 

(22b) 


energy is excited for^ithe/He^^^ that no elec tronic 

1000 K) and that the rotational starf is a 8 ° od assumption below 

vibrational modes are at the eauil^ri. "* #re fully excited whil « the 
Since the characteristic vibrational t-!!! exc * tatlon for a harmonic oscillator, 
the latter assumption is l°* * 2 U ^ (6 v » ^159 K) 

rotational temperature for H 2 is rather larc« h ° W o= e fl fc J® cha raccer istic 
motion of equation (22b) can be expected to 8 hold r " 85 ’ 42 K) * 80 the a *>P r °xi- 
temperature end of the range near 100 k 7 * M 8ome error at the low ~ 

The thermal conductivity of the gas mixture 1, now expreeaed 


5n 

32 
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to^he^eanbfree t path^ln h the^pure^lth D ' ean P<,tl ' f ° r ener *V «->■>■>- 

sta? sl-sss^ i f :r ^ 

about equal, which is the case for the^res^r^ ’ lf ^ molecular masses are 
tion of excess kinetic energy! assoc T** ° f lnterest ‘ The frac- 
centers at impact, transferred between 
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ik 


4m /m, 
i k 

1 + m /m, 
i k 


(24) 


Accordingly, the excess kinetic 
uct of the collision rate 


energy transferred 
and the fraction 


per unit time is the prod- 
fik summed over all k 



( 25 ) 
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and the ratio 




can be expressed 


X ° 


2m V-i 
*i Q u P * "A 

k*i 2 
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(26) 


in this approximation, the thermal conductivity of the mixture is given by a 
itie^o^mixture^ 68810 ” analogoU8 to the Unear *zed formula for the viscos- 


X k * (»*•) 


(27) 


Undoubtedly, an empirical correction term, such as that provided by 
“ d " llke f °' the viscosity of fixtures, could accountfo L 
difference in the cross sections for diffusion and for energy transfer and 
thus improve the theoretical moeel somewhat; however, such a correction has 
not been established st present. Since the molecule^ masses ace wany aqua! 

iercmt af 8 K’ f tdtfectlon could only be a modest one, of the order of a few 

percent at best. Such a correction takes on major importance only in ionized 

energy1rlns£r.'’ e,> ' ll8htVel8ht eUc “°“ “ r ry a large fraction* t£ 


Thermal Conductivity of Pure Molecular Hydrogen 

aK1 Th ® coefficient of thermal conductivity of pure H 2 is reproduced reason- 
ably well over the temperature range of 100 to 1000 K bv the formula 

0 /2T 

V 


k = /2m^y/‘fl9 + / °y /2T \ ^| R 

12 \ * / y slnh 0J2l) I mQ 


1 . 10 X 10 -- <&.&>> !' 6m3 


( 1 + 


'300)0 .6983 H, / 3Q8Q/T \ a| 

*P) L 4 U^)j 


cal /cm deg sec (28) 


where again m is the mass of the H; molecule. Accordingly, the energy 
transfer cross section Q, for Hj - H 2 collisions, is taken to be 


Q - 2«.l»10-*« (M)°- »»•*(, + 


(29) 


The deviation of equation (28) from the H ; > coefficients of thermal con- 

devJnM ty *' Vetl b y« 11 f enr «th et al. (ref. 6) is shown in figure 3. The 
deviations, generally less than 1%, are adequate for present purposes. No 
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Figure 3.- Deviations of equation (28) from published thermal conductivities 

of pure Ho 


doubt, some of the deviation near 100 K is due 
the equilibrium rotational excitation. 


to the approximate treatment of 


Ihermal Conductivity of Pure Helium 


2 ^2mkT y /? 


15 R_ 
4 mQ 


3.42.10-^ (T/3q0)^- 74 


(1 - 13.7/T) cal/cm deg sec (30) 


cross section ‘q, ' f o^^^e'coiHsi^; if tokenT-V''' em ’ rSy tra " 8f '' r 


5.96*1O- 16 0~\ 


0 . 1 ’ 4 1 ; 


(- - l ¥)~' 


til) 


dm- 1 J'u'vfi vf tMirr.,:;! 0 ! , ::r. <* t»«i ««. 

deviations are senerallv within' 27 I ?! ore “ hown ln fl 8":e 4. The 

hy Powell «ref.\o> are’ aTf'sh owf wISTjS “ tlT, 

mixture formulas used for therma! coXi^ni!!!! ‘"--rtainties «f the 
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o POWELL. HO, LILEY 
(REF. 9) 

9 POWELL (R6F. 10) 


0 200 400 W m ^00 

TEMPERATURE, K 

Figure «.- Deviations of equation (30) from published themal conductlvitiee 

of pure He 

Thermal Conductivity of Jovian Atmospheres 

a function of tLperftu^e V in y taMe°3 e L^°the n t atm08Phere mixtures ls shown as 
1000 K. Again, the first col»™ < \t \ temperature range from 100 to 
thermal conductivity of pure Ho (which 6 emperat “ re 5 the second column is the 

Model and the Warm Light^ode^of^ tte J^T* 8 u° T" the C ° o1 Li * ht 
the conductivity of an 89% H u 4 tmosphere) , the third column is 

Nominal Jovian It^ospSere Model) f nfthe f^! ^ C ° rre8ponds to the 
a 78% H 2 - 22% He mixture (wMci/,,* ° th j inal column is the conductivity of 

Warm Heavy Model of the Jovian Atmosoh^e^ 8 Th both . the C ° o1 Heav y Model and 
Model Atmosphere (with 11% He) is fl iL *■ ® on< * uc tivity of the Nominal 

independent of tenperature. for pure H,, 

sphere Models (vith 222 He) fa „t , y * t * le con ductivity of the Heavy Atmo- 
ivith 222 He) is almost constantly A2 less than for pure H 2 . 
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TABLE 3.- THERMAL CONDUCTIVITY OF JUPITER MODEL ATMOSPHERES 


i 


Thermal conductivity, mcal/cm-deg-sec 


Temperature, K 

Light atm. a 


100 

0.1623 

150 

.2422 

200 

.3156 

250 

.3841 

300 

.4486 

350 

.5099 

400 

.5685 

450 

.6249 

500 

.6794 

550 

.7323 

600 

, .7839 

• .8344 

.884 

650 

700 

750 

.9331 

800 

.9816 

850 

1.0299 

900 

1.0779 

950 

1.1258 

1000 

1.1737 


Nominal atm.^ 

Heavy atm/ 

0.161 

0.1597 

.2395 

.2365 

.3115 

.307 

.3787 

.3727 

.4419 

.4345 

.5021 

.4933 

.5596 

.5496 

.6149 

.6038 

.6684 

.6561 

.7203 

.7069 

.771 

.7565 

.8205 

.805 

.8693 

.8527 

.9174 

.8998 

.965 

.9464 


1.0123 

1.0594 

1.1064 

1.1533 


.9926 

1.0387 

1.0845 

1.1303 


‘*Thennal conductivity of pure H 2 . 

^Thermal conductivity of 89% H 2 - 11% He mixture. 
1 Thermal conductivity of 78% Ho - 28% He mixture. 


Table 4 shows the thermal conductivity of the Jupiter Nominal Atmosphere 
as a function of altitude and the ratio of the mixtuce conductivity to the pure 
Ho conductivity. 
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TABLE 


TEERMAL CONDUCTIVITY OP JUPITER NOMINAL ATMOSPHERE 


Altitude, km 


- 270.79 

- 246.93 

- 224.62 

- 203.74 

- 184.23 

- 165.98 

- 148.91 

- 132.96 

- 118.03 

- 104.09 

- 91.06 

- 78.87 

- 67.49 

- 56.86 

- 46.94 

- 37.66 

- 29.03 

- 20.97 

- 13.48 

- 6.49 

0 

6.02 

11.6 

16.75 

21.43 

25.79 

29.97 

34.04 

38.09 
42.16 

46.32 

50.57 

54.91 

59.36 

63.9 


Temperature, K 


681.77 
637.33 
595.74 
556.83 
520.44 

486.39 

454.55 

424.78 
396.93 
370.89 

346.51 

323.66 

302.25 

282.17 

263.31 

245.59 

228.92 

213.21 

198.4 

184.43 

171.23 

158.77 

147 

134.3 

121.6 

116.65 
111.7 
110.9 
110.1 
112.35 

114.6 

117.29 

119.98 

122.66 
. 125.35 


Conductivity, 
mcal/cm— deg sec 


0.8516 

.8081 

.7667 

.7273 

.6898 

.654 

.6199 

.5873 

.5561 

.5264 

.498 

.4707 

.4447 

.4198 

.3959 

.3729 

.3509 

.3297 

.3093 

.2897 

.2707 

.2525 

.2349 

.2156 

.1958 

.188 

.18 

.1788 

.1775 

.1811 

.1847 

.189 

.1932 

.1975 

.2017 


Conductivity 
ratio, mix/H 2 


0.9633 

.9834 

.9835 

.9836 

.9837 

.9839 

.984 

.9842 

.9843 

.9845 

.9847 

.9849 

.9851 

.9854 

.9856 

.9859 

.9862 

.9866 

.987 

.9874 

.9878 

.9883 

.9888 

.9895 

.9903 

.9906 

.991 

.991 

.9911 

.9909 

.9908 

.9906 

.9904 

.9902 

.99 


TABLE 4.- CONTINUED 


r 

. 1 1 

Altitude, km 

Temp era ture , K 

Conductivity, 

Conductivity 




mcal/cm-deg sec 

ratio, mix/H 2 


128.04 

130.73 

133.41 


0.2059 

.2101 

.2142 


0.9898 

.9897 

Qft<K 


93.2 

140.09 

98.37 

142.09 

103.62 

144.09 

108.94 

146.08 

114.34 

148.08 

119.81 

150.08 

125.36 

152.08 

130.98 

154.07 

136.67 

156.07 

142.44 

158.07 

148.28 

160.06 

154.2 

162.06 

160.19 

164.06 

166.25 

166.05 

172.39 

168.05 

178.61 

170.05 

184.9 

172.04 

191.26 

174.04 

197.7 

176.04 

204.21 

178.03 

210.8 

180.03 

217.46 

182.03 

224.2 

184.03 

231.01 

186.02 

237.89 

188.02 

244.85 

190.02 

251.89 

192.01 

259 

194101 

266.19 

196.01 

273.45 

198 





TABLE 4.- CONTINUED 


Altitude, km 


Temperature, K 


Conductivity, Conductivity 

mcal/cm-deg sec ratio, mlx/H 2 


280. 78 
288.16 
295.53 
302.91 
310.29 

317.67 

325.05 
332.44 
339.82 
347.21 

354.6 
362.39 
370.98 
380.37 
390.57 


200 

200 

200 

200 

200 

200 

200 

200 

200 

200 

200 

221.67 
243.33 
265 

286.67 


0.3115 0.9869 

.3115 .9869 

.3115 .9869 

.3115 .9869 

•3115 .9869 


.3115 

.3115 

.3115 

.3115 

.3115 


.9869 

.9869 

.9869 

.9869 

.9869 


•3115 .9869 
•341.2 .9864 
• 37 .986 

• 398 .9856 
.4254 .9853 


401.58 

413.39 

426 

439.42 

453.65 


308.33 
330 

351.67 

373.33 
395 


.4522 

.4784 

.504 

.5292 

.5539 


.9851 

.9848 

.9846 

.9845 

.9843 


468.69 

484.54 

501.19 

518.66 

536.95 


416.67 

438.33 
460 

481.67 

503.33 


.5783 

.6022 

.6258 

.649 

.6719 


.9842 

.9841 

.984 

.9839 

.9838 


556.04 

575.96 

596.69 

618.24 

640.6 


525 

546.67 
568.33 
590 

611.67 


.6946 

.7169 

.739 

.7609 

.7826 


.9837 

.9837 

.9836 

.9835 

.9835 


.4 


663.79 
687.8 
712.64 
738.3 

764.79 


633.33 
655 

676.67 

698.33 
. 720 


.8041 

.8255 

.8466 

.8677 

.8886 


.9834 

.9834 

.9833 

.9833 

.9832 
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CONCLUDING REMARKS 


The approximations used in estimating the viscosity of Jovian atmosphere 
mixtures are believed reliable within a few percent. The approximations used 
in estimating the thermal conductivities of Jovian atmosphere mixtures hold 
uncertainties of the order of 10 to 15%. In both cases , the mixture trans- 
port coefficients are almost a constant factor times the value for pure 
independent of temperature. This permits a quick, easy method of evaluating 
the appropriate values for the mixtures. 

The transport coefficients included here are appropriate only for sub- 
sonic flow fields in the Jovian atmosphere. At high supersonic speeds, the % 
dissociates and both the He and H atoms are electronically excited and ionized 
in flow following strong shock waves and in boundary-layer flow. In this 
case, the dissociation and ionization reactions promote very large peaks in 
the thermal conductivity over the range of pressures and temperatures where 
the reaction rates are significant. The ionization of the species also 
decreases viscosity by an order of magnitude or more, because the long-range 
coulomb forces prevent long mean-free-path momentum transfer in the gas. The 
characteristic dissociation temperature of H 2 (that is, the ratio of dissocia- 
tion energy U to the Boltzmann constant k) is 51,965 K, and the characteristic 
electronic excitation and ionization temperatures are several times larger. 
These reactions will not become Important until temperatures are of the order 
of 0.1 the characteristic temperature. In other words, the mixture formulas 
used here could be applied to temperatures of several thousand degrees. 

However, in such a case, the simple formulas used for the viscosities and 
thermal conductivities of the pure H 2 and He would have to be modified; they 
were designed only to apply in the range from 100 to 1000 K. 

Tables 5 and 6 are printouts of the BASIC language programs that have 
been used to calculate the coefficients of viscosity and thermal conductivity 
of H 2 “ He mixtures. In the future, we may well find that other mixtures are 
appropriate for the Jovian atmosphere, in which case the appropriate values of 
the mole fractions xj (of H 2 ) and X 2 (of He) can be introduced as input to 





TABLE 5- BASIC PROCRAM POR VISCOSITY OP JOVIAN MODEL ATMOSPHERES 


VISCOSITY OF JUPITER MODEL ATMOSPHERES, MICROPOISE' 


10 PRINT "TABLE I - 
20 PRINT 

40 PRINT " T ’ DEG K "’" LIGHT A™-, "NOMINAL ATM", "HEAVY ATM" 

! J! 0 !' ^ raction of H2 in nominal atmosphere model 

i m°i* f raction of He * nominal atmosphere model 

t m°J 5 raCt j° n ° f H2 ** 1 heay y atmosphere model 

f TnJ*f r ? Cti ° n ° f He ln heavy atmosphere model 

i initial temperature, deg K 


! 


t 


50 XI". 89 
60 X2". 11 
70 Yl". 78 
80 Y2-.22 
90 T"1E2 

100 FOR I"1 TO 19 
140 Nl»90 . 6* (T/300) A . 6658 
142 N1»N1/ (1+4/T) 

150 N2-191. 6* (T/300) A . 7176 
152 N2«N2/ (1-11. 4/T) 

160 Ql*32 . 3* (1+4/T) * (300/T) A . 1658 

H2 < * 2 ° 21 • 5* (1-11 . 4/T) * ( 300/T) A . 21 76 
180 Q3»(Ql A .5+Q2 A .5)/2 
190 Q3"Q3*Q3 

200 L1-1+.7967*(X2/X1)*(Q3/Q1) ! 

220 N-N1/L1+N2/L2 1/X ? > VUm Q ^ f ' f L <P“ r e)/L(,»ix) for He 

230 L3«l+.7967*(Y2/Y1)*(Q3/Qi) y ° ^Ratio ^? 08 ^?/? '< ® iaro P olse 

iss Potto sm s 2 

280 NEXT I NeXt value of temperature, deg R 

999 END 


Viscosity of H2, micropoise 

Viscosity of He, micropoise 
f Cross section H2-H2, Ang2 
I Cross section He-He, Ang2 

! Cross section H2-He, Ang2 
Ratio L (pure) /L (mix) for H2 
Ratio L(pure)/L(mix) for He 

«•( «pi1 » 1 * 
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TABLE 6.- BASIC PROGRAM FOR CONDUCTIVITY OF JOVIAN MODEL ATMOSPHERES 


10 PRINT "TABLE III - THERMAL CONDUCTIVITY OF JUPITER MODEL ATMOSPHERES" 
U PRINT " MILLICAL/CM-DEG-SEC" 

20 PRINT 

30 PRINT "T,DEG K", "LIGHT ATM" , "NOMINAL ATM", "HEAVY ATM" 

40 PRINT 

! Mol fraction of H2 in nominal atmosphere model 

! Mol fraction of He in nominal atmosphere model 

! Mol fraction of H2 in heavy atmosphere model 

! Mol fraction of He in heavy atmosphere model 

! Initial temperature, deg K 
100 FOR 1-1 TO 19 

140 K1-. ll*(T/300) A . 6983/ (1+49. 4/T) 

142 V-3079.5/T 

144 V»2*V/(EXP(V)-EXP(-V)) 

146 K.1-K1* (4 . 75+V*V) ! Conductivity of H2, millical/cm-deg-sec 


50 XI- . 89 
60 X2-.ll 
70 Y1-.78 
80 Y2-.22 
90 T-100 


! Conductivity of He ,mcal/cm-deg-sec 
! Cross section H2-H2, Ang2 
J Cross section He-He, Ang2 


150 K2- . 3418* (T/ 300) A . 7412/ (1-13. 74/T) 

160 Ql-26. 1* (1+49. 4/T) * (300/T) A . 1983 
170 Q2-5. 96*(1-13. 74/T)*(300/T) A . 2412 
180 Q3- (Q1 A . 5+Q2 A . 5) /2 

IaS ??“? 3 * Q3 1 Cros ® section H2-He, Ang2 

200 L1-1+. 7698* (X2/X1)*(Q3/Q1) ! Ratio L (pure) /L (mix) for H2 

oon ^“| , ‘t 1 -* 088 ^*^ x i/ x 2)*(Q 3 /Q2) ! Ratio I. (pure) /L (mix) for He 

'ZzZ K “ K 1/L1+K2/L2 ! Thermal conductivity nominal atm,mcal/cm-deg-sec 

230 L3*l+. 7698* (Y2/Y1) * (Q3/Q1) | lutfo L(pure)/I.(mlx) for B2 

240 L4-1+1 . 0887* (Y1/Y2) * (Q3/ Q2) ! Ratio L (pure) /L (mix) for He 

bbtS^ 3 ^ 2 /^ 4 * Thermal conductivity heavy atm, mcal/cm-deg-sec 

260 PRINT T,1E-4*INT(1E4*K1),1E-4*INT(1E4*K),1E-4*INT(1E4*K3) 

270 T-T+50 I Next value of temperature, deg K 
280 NEXT I * 

999 END 
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